New method for studying steady states in quantum impurity problems: The 

interacting resonant level model 
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We develop a new perturbative method for studying any steady states of quantum impurities, in 
or out of equilibrium. We show that steady-state averages are completely fixed by basic properties 
of the steady-state (Hershfield's) density matrix along with dynamical "impurity conditions". This 
gives the full perturbative expansion without Feynman diagrams (matrix products instead are used) , 
and "re-sums" into an equilibrium average that may lend itself to numerical procedures. We calculate 
the universal current in the interacting resonant level model (IRLM) at finite bias V to first order 
in Coulomb repulsion U for all V and temperatures. We find that the bias, like the temperature, 
cuts off low-energy processes. In the IRLM, this implies a power-law decay of the current at large 
V (also recently observed by Boulat and Saleur at some finite value of U). 

PACS numbers: 73.63.Kv, 72.15.Qm, 72.10.Fk 



Impurity models describe mesoscopic quantum objects 
in contact with large conducting leads. Quantum dots 
are much-studied examples, and experiments in a bias 
voltage [H lead to accurate descriptions of their non- 
equilibrium steady-state properties. In impurity models, 
such states are of high interest as they pose the theoreti- 
cal challenge of capturing the effect of non-equilibrium in 
a truly quantum system, yet they constitute the simplest 
non-equilibrium situation, where properties are time in- 
dependent. At low energies, interactions between the 
leads' Landau quasi-particles and the impurity occur in 
the s-wave channel, and the spectrum may be linearised 
around the two Fermi points. Then, the universal be- 
haviour is described by free massless fermions on the half 
line (bulk conformal field theory) with a non-conformal 
boundary interaction at the end-point. Many methods 
are known for studying equilibrium behaviors in such 
models, but understanding and accessing properties of 
non-equilibrium steady states is still a much harder task. 
Wilson's picture does not apply, and, for instance, the ef- 
fect of a bias on low-energy processes (which determines 
the large-bias current) is still under study 0, Q • 

Two calculational schemes exist: the real-time 
Schwinger-Keldysh formulation and the scattering-state 
Lippman-Schwinger formulation. The former relies on 
the infinite extent of the half line in order to absorb 
the energy released by relaxation from an appropriate 
"uninteracting" density matrix to the steady state. It 
requires relaxation mechanisms (see, e.g. [3[), whose ab- 
sence may lead to pathologies in perturbative expansions 
The latter describes directly the expected end result 
just from "how the state looks" asymptotically far from 
the impurity. Hershfield's Y operator Q (see also the 
studies 0, & 0) gi yes a "steady-state density matrix" 
that encodes these scattering states. This is interesting, 
since a non-equilibrium steady state is not described by 
the usual density matrix, but it is still hard to apply to 
interacting systems. Exact results for integable models 



occur when the exact quasi-particles do not mix the bi- 
ased particle baths [§]. A recent proposal || suggested, 
in the interacting resonant level model (IRLM), a free- 
dom in choosing the quasi-particles to get around this 
restriction, but the construction of Bethe ansatz eigen- 
states still needs justifications (for another recent study 
of the IRLM, see 0). 

In this paper, we develop a new method for study- 
ing steady states in quantum impurity models from ba- 
sic properties of Hershfield's Y operator and "impurity 
conditions" . Impurity conditions are part of the equa- 
tions of motion that relate local operators on both sides 
of the impurity. In integrable models, they fix the scat- 
tering matrix of Bethe ansatz, and they are used in 11 1 
to describe the infrared fixed point of the IRLM. We use 
them to develop an efficient perturbation theory in the 
scattering-state formulation. They allow us to bring chi- 
ral fields (representing collective massless modes in the 
bulk) through the impurity to the side where Hershfield's 
steady-state density matrix looks like the usual equi- 
librium one. The steady-state condition is then solved 
iteratively in terms of equilibrium averages, giving the 
full perturbative expansion without Feynman diagrams, 
where the combinatorics is dealt with using matrices on 
the space of impurity operators. 

This is an improvement over Keldysh perturbation the- 
ory - there is no special time ordering; and over other ap- 
proaches in the Lippman-Schwinger formulation - there 
is no need for explicit scattering states or Y opera- 
tor. It gives full perturbative expansions with respect 
to marginal or marginally relevant operators, which is 
impossible by standard methods. The re-summed ex- 
pansion is an equilibrium average that may lend itself 
to numerical techniques. The method is also related to 
"equations of motion methods" [13] ; physically motivated 
truncations give non-perturbative results. Conceptually, 
the method shows how basic properties of the steady- 
state density matrix are sufficient to uniquely determine 
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quantum averages; this may help in understanding vari- 
ational approaches out of equilibrium. 

In the IRLM, we calculate the universal non- 
equilibrium current for all voltages and temperatures to 
first order in Coulomb repulsion U > (we verify that 
Callan-Symanzik equations hold). We find that large dif- 
ferences between the leads' filling energies (determined by 
the voltage) and the impurity energy cut off low-energy 
processes, giving a power-law decay at large voltages. 
Such a decay was also observed recently at a particu- 
lar non-zero value of U [H}. We note that the Bethe 
ansatz eigenstates of @ do not satisfy our impurity con- 
ditions, except at U = 0, oo (where exact quasi-particles 
are free). Finally, we give an argument to obtain the 
current at U = oo, and find agreement with Q and with 
free-theory expectations. 

The hamiltonian. We use units where e = h = 1. 
The IRLM is described by a quantum field theory hamil- 
tonian H = Hq + Hj. Hq represents two baths of one- 
dimensional free massless, spinless fermions (with quan- 
tum fields tpi(x) and if>2(%)) on the half line, and Hj the 
interaction with the impurity, an energy level that can be 
singly-occupied, at the end of the half line. Since there 
is no back-scattering, we can "unfold" the two baths and 
get two chiral right-moving fermions on the line. Using 
V^tpe — V^i + V^ipo — ipi — 4>2, the hamiltonian is 



i / dx (Tpl{x)dMx)+ipl{x)d x i> (x)) (1) 



Hi = t(iPl(0)d+dHe{0)) + U(n e (0)+n o (0))D + e d D 

where D = d)d and n eM (x) = t/>| (x)ip e>0 (x). Op- 
erators satisfy canonical anti-commutation relations 
{^ e {x)^t(x')} = {ip o (x),ipt(x')} = S(x - x 1 ) and 
{d,d^} = 1, d 2 = (d^) 2 — 0. U represents the Coulomb 
repulsion between the leads' fermions and an occupied 
impurity, is the impurity energy level, and t is the hy- 
bridization, which essentially determines the number of 
current-carrying states available. It is useful to keep in 
mind that in the unfolded picture, x < means before 
interacting with the impurity, and x > means after. 

Steady states. In quantum impurity systems, steady 
states are also time-independent quantum states. They 
arc defined by prescribing the baths to be, at x < 0, in 
thermal equilibrium at temperature T with a chemical 
potential V associated to a local conserved charge Q of 
.Ho, and by asking for time- independence with respect to 
the dynamics given by the full hamiltonian H. That is, 
quantum averages (• • •) = Tr (p- ■ •) /Tr (p) are described 
by the density matrix [5| p = exp [—(H — VY)/T], where 
Y looks like Q before interacting with the impurity: 



x 1 <0,x 2 <0,. 



(0 1 (x 1 )0 2 (x 2 ) ■ ■ ■) 
Tr (exp [-i(ff - VQ)] 1 (x 1 )0 2 (x 2 ) 



(2) 



Tr (exp [-±(H - VQ)]) 
(Oi(x) are local operators), and is conserved by the dy- 



namics 17]: [H, Y] = 0. As we will see, these two 
properties are sufficient to uniquely characterise steady- 
state averages. For our purposes, the charge Q = 
(N 2 - Nx)/2 with Ni = f dxipl(x)ijji(x) and the current 
J = — i[H, Q] — tRe(id^tp o (0j) describe a bias voltage 
and the associated particle current. 

The method. Impurity conditions. Eigenstates of H 
can be constructed by linear combinations of states of the 
form Jdxdx'--- g{x,x',...)iPl {x)^l {x r )--- (<?) 0,1 |0) 
with pseudo-vacuum |0), ip e<o (x)\0) — 0, d\0) = (a 
vacuum for both the fermionic Fock space T and the im- 
purity space I) . The Hilbert space Ti comprises states up 
to an energy distance A from the pseudo-vacuum - this 
represents finite leads' bandwidths of order A, a part of 
our cutoff scheme. Because of the linear spectrum in 
the bulk, the first-quantised wave functions g have fi- 
nite jumps at the impurity site. Hence, (v\ ty(x) \w) has 
a jump at x = 0, where \v), \w) € Ti and ^{x) is an 
ultra-local field (no derivatives) at x times impurity op- 
erators. Using [H ,^(x)] = id x ^{x) in (v\ [H,Hf(x)) \w), 
the jump is, as an operator on Tt (with a < 0, b > 0): 




[ dx - [ dx] [#j,$(a 

Ja J0+ J 



(3) 



Eq. ([3]) has ambiguities: Hj contains operators at 
x = 0, where there are jumps. This is an artefact of 
the un-physical point-like nature of the impurity in the 
model, and is lifted by "spreading" it in a region around 
x = (another part of the cutoff scheme). Equivalcntly, 
one can "resolve" the hamiltonian, by only using lim- 
iting operators ^ e ,o(0 + ), ipe,o{0~)- A natural symmet- 
ric resolving [l8j|. which we will use, is the replacements 
^,0(0) !— ► (i/ , e ,o(0 + ) + ip e , o (0~))/2 (and hermitian con- 
jugate) in Hj. In a longer version of this work we 
will discuss these issues and construct the Bethe ansatz 
eigenstates of the resolved hamiltonian. 

Taking, then, ty(x) = ip e (x) in (|3|) gives 



1 + ^D),-,((P 



-itd 



and a similar equation holds with tp e 1— > ip Q and t 1— > 0. 
From ((3J, these hold only as operators on Ti. However, it 
can be explicitly checked on Bethe ansatz eigenstates that 
they hold more generally as linear maps H. — > T®X. This 
is the meaning of the equality symbol above. This fact 
may be due to integrability of H , and fixes the scattering 
matrix of Bethe ansatz eigenstates (different from that of 
[I if U ^ 0,oo). It also gives "right-limits" in terms of 
"left-limits" . For instance, with u = 



2iU 
2i-U ' 



ipe(0 + ) = -M + (1 - iuD) ip e (Q-) 



(4) 



Other equations are obtained from ((4]) by pre- 
multiplication by d, d< and ipl(0 + ), ipl(0~) and by her- 
mitian conjugation and replacement e 1— ► o, 1 1— > 0. These 
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are our impurity conditions; they stay valid when multi- 
plied by local operators at i ^ 0. What they tell us is 
how local operators that are just to the right of the im- 
purity can be brought to its left only using the dynamics. 
They give for the current J = iRje(i(Fip o (0~)) on H. 

Steady state conditions. Using impurity conditions, the 
steady state condition (\H,ty(x))} = for x < can be 
written with all local operators to the left of the impurity. 
We write ^(x) — djO(x) for some j = 1,2,3 where we 
denote b\ = d, 62 = d \ 63 = D. Then, we have 

~id x {b 3 G(x)) = l A 3 (b,O(x)) + ((c 3 +Y,b t E U3 (0-))O(x)} 



through bosonisation techniques and give a perturbative 
expansion in t. 

Divergencies. The expansion ([7]) is plagued by linear 
and logarithmic divergencies as A — > 00 . Linear divergen- 
cies (from normal-ordering with respect to Hq in Eij(x), 
n(x) 1 — y. n(x) :, and from collisions : n(x) : : n(x') : as 
x ~ x' in 10) will be absorbed into a re-definition of 
e d [l4[. Logarithmic divergencies are accounted for by 
the renormalisation group - a good alternative regulari- 
sation that gives them is simply to "avoid" the impurity 
by changing the upper integration limits in to — A -1 . 

Results. We take U > 0, and that regularisation. We 
find the current (divergent and finite parts as A — * 00) 



where Aj is %^ + ie d , -y — ie<j, t 2 and Cj is — ttjj e (0 ), 
tipl(Q~), for j = 1,2,3, and 



(J) =m[I cs + U 



I r m 



d_ 

dm 



Ir 



D 



0{U 2 ) 



-un(x) — ttpl(x) 
Ei_j{x) = I un{x) —t"ip e (x) 

ituip e (x) ituipl(x) 



(5) 



1,3 



with n = n e + n a . This is valid also with 0(x) 1— > 0(x + 
xi) ■ ■ ■ 0(x + x n ) for x + xj < 0. 

Full perturbative expansion in U. This differential 
equation can be integrated, choosing integration con- 
stants for a finite limit x — > —00 of (bjO(x)): 

(bjO{x)) = ie- A * x f dx'e A ' x ' (( Cj +J2 b iE id (0-))O(x')) 

(6) 

Solving ([6]) by iteration, with J^i biE iy j (0 )0(x') in place 
of bjO(x) at x — > 0~ and repeating, gives: 



3 
i=i 



dxo ■ ■ ■ dx n {c T e n E{x n ) x 



l E{ Xl + ... + x n )e Ax °O(x + ...+ x n ))° 

(7) 

where A — diag(Ai, A2, A3), c T = (01,02,03), E is the 
matrix (0, and 6 = (d, 2 , 3 ) T . All terms of this 
series are averages of local operators strictly to the left of 
the impurity. Hence (0 was used, and (• • -)° means free 
equilibrium averages like on its right-hand side. This is 
our main result: we have obtained the full perturbative 
expansion. Essentially, the use of impurity conditions 
has simplified the combinatorics involved. For [7 = 0, 
the first interation gives the current in closed form, in 
agreement with e.g. @. 

The expansion (0 can be written using path- 
ordered integrals as i f_ dx (c T- Pexp f£ dx' (-iE(x') + 
A)0(x))°. This can be evaluated using appro- 
priate random variables instead of operators, which 
may be useful for numerics. It can also be 
re-summed using Baker-Campbell-Hausdorff formula, 
iJ° 00 dx(c T e- t( - Ho+lA+E ^d(0)e lHoX )°, where E = E(0). 
This should allow a full account of the diagonal part of E 



where m = t 2 /2, I cc = \ £ ±> ±, F(0, m ± ie d ±' iV/2) 
L(A), I c = ±£ ± F(0,2m'± iV/2) + L(A), I cs 
i J2 ± F(m ± ie d + iV/2, m±ie d - iV/2) and 



B = mT' 



dy 



' , e m ( y+z ^cose d (y- z) , 
dz — — ^ -(y - z)X 



Vz ■ V(y—z) 

cos sin — ^ — '- 
sinh7rT(y — z) 



sinh irTz 

Vz Vy 



Vy 



y cos sin + z cos sin X^- 
(y — z) smhuTy 



with L(A) = ^log|^ and, with ip(z) = dla.T(z)/dz, 
F( 0|6 ) = I(^(_^+i)-^(_^+i)). 

The average current satisfies the Callan-Symanzik 
equation (^m^ - A^) (J) = 0. This checks RG out 
of equilibrium to all orders in m. It means that the scal- 
ing limit where the bandwidth A is greater than all other 

scales is taken with the parameter T K " = mA * fixed 
(to first order in U). This trades the microscopic cou- 
pling m for the scale Tk that characterises the scaling 
limit. The renormali^ed result is as above with the re- 
placements m 1— » T K " and Ai->/i, for some "per- 
turbative scale" n. Since the parameter U does not flow 
[Tot HBj , the renormalised perturbative expansion is au- 
tomatically valid for all T, V, e d , Tk -C A, i.e. in the uni- 
versal regime. Note that the exponent 1 + U/n + 0(U 2 ) 
is in agreement with standard bosonisation results. 

The full result is independent of fj,, which is only cho- 
sen to keep perturbative coefficients finite. Then, fi is 
the greater of Tk or any physical infrared (IR) cut- 
off (quantities that do not allow arbitrarily many low- 
energy processes, like the temperature). It turns out 
that with (j, - V( A + + T + t k){&- + T + T K ), with 
A± = \V/2 ± e d \, when any of T,T#,A-t is/are much 
larger, the first order perturbative coefficient stays finite. 

Why are the energy separations A± cutting off low- 
energy processes? Recall that transfers of particles occur 
around a point in space, so that by Heisenberg's uncer- 
tainty principle all energy scales in the leads may be in- 
volved. When at least one of A± is high, there cannot 
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FIG. 1: Universal current vs. voltage, Tk fixed, T = = 0. 
A universal feature is the large- V decay with a [/-dependent 
power-law. At small V, its decrease with U is not universal, 
and depends on the definition of Tk ■ 

be many low-energy transfers without many high-energy 
transfers on at least one side of the impurity. Here we 
find that each side contributes a factor of a square-root 
to the IR cutoff. That the voltage is a good infrared cut- 
off was also observed perturbatively in the Anderson [l6| 
and Kondo 0] models. 

One implication is that the current vanishes as a 
power law as F » T,T K ,\e d \ (see Fig. 1). At 
small temperatures, V 3> Tk ^ T, |e<j|, we have 

(J) ~ \T K (8eW/ 2 ^)" + ° {U2) (1 + 0(U 2 )). This is 
explained as follows. Observe that at U = 0, the current 
saturates at large V, to a value proportional to the num- 
ber of current-carrying hybridized states, roughly m/A. 
According to RG, the effect of U ^ is that at large 
energies, m/A decreases with a power law. Since V is an 
infrared cutoff, at large V the system is at large energies. 
Note that if | | is kept of the same order as V/2, the 
power law changes to . 

On the other hand, the low-T linear current, Tk ^> 
T, V, \e d \, is (with f = ttT/Tk) 2n(J)/V ~ 1 + g 2 f 2 + 

g 4 T 4 + 0(T 6 ) with g 2 = - (l&e^/Vy * /3 and univer- 
sal ratio g^/g 2 = 21/5 — U/n, everything up to 0{U 2 ). 
The universal ratio is in agreement with IR conformal 
perturbation theory [ill ]. 

Note that fixing Tk, the current decreases as U is in- 
creased, for all voltages and temperatures (see Fig. 1 for 
the zero-temperature case). The fact that Tk can be cho- 
sen in such a way is universal. The opposite is observed 
for finite but large enough bandwidth, finite voltage and 
fixed t, in agreement with (il. Ho|. 

Physical arguments also give the current at U — oo, 
where u = 2i. Here, there cannot be a fermion at x = 
and an occupied impurity level simultaneously. Hence, 
both operators dn(0~), d^n(0~) must average to zero. 
In ([6]) (with linear divergencies substracted) , this brings 
the matrix E to its U — form, so that the current is 
m I cs , in agreement with Q . 

Conclusions. We have developed an efficient method 
for studying steady states in impurity models, applied it 



to the current in the IRLM, and identified the non-trivial 
voltage-related quantities (A±) that cut off low-energy 
processes. It can be applied to any single-impurity model 
with linear spectrum and with left-/right-moving sepa- 
ration of the leads' quasi-particles (where unfolding can 
be done): Kondo and Anderson models, and impurity 
models where leads are Luttinger liquids (bosonic col- 
lective modes separate). In the non-equilibrium Kondo 
model with an impurity magnetic field, it may simplify 
the treatment of the pathologic perturbative expansion 
[3] • Although our method may be restricted to integrable 
dynamics (this is to be clarified), it does not rely on inte- 
grability of the steady state; the non-equilibrium steady 
states in the Kondo and Anderson models may not be 
integrable. It would be very interesting to develop ef- 
ficient, out-of-equilibrium numerical methods from the 
re-summed perturbative expansion obtained here. 
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